The scalar and tensor glueballs in the valence approximation* 
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We evaluate the infinite volume, continuum limit of and 2 glueball masses in the valence approximation. 
We find mo++ = 1740 ± 71 MeV and m^++ = 2359 ± 128 MeV, consistent with the interpretation of /o(1710) as 
the lightest scalar glueball. 



Encouraging results obtained recently for light 
hadron masses in the valence approximation [1] 
have prompted us to investigate the glueball spec- 
trum. For both 0"*""*" and 2+"'" states, a quark- 
antiquark system has odd orbital angular momen- 
tum. The resulting angular momentum barrier 
should tend to suppress the quark-antiquark an- 
nihilation required for mixing with glue states. 
An argument [2] based on the observed near de- 
generacy of corresponding pairs of isovector and 
isoscalar mesons composed of a quark and an an- 
tiquark with nonzero orbital angular momentum 
supports this expectation. Thus we believe there 
is a reasonable chance that our results provide 
fairly reliable predictions for the real world. 

We use square loop glueball operators con- 
structed from smeared links. Smearing eliminates 
some of the high frequency noise and increases the 
projection to the ground state [3]. Although dif- 
ferent methods have been proposed and used in 
the past, most of them consist of defining smeared 
5{7(3) link matrices by adding and multiplying 
together neighboring links in a gauge invariant 
way [4-7]. Part of our work uses the operators 
proposed in Ref. [7], in which the smearing of 
Ref. [4] was used to construct arbitrarily large 
Nl X Nl loops. In Ref. [7] it was also found 
that this method is similar in performance to 
the smearing proposed by Tepper [5] and by De- 
Grand [6]. Our operators are then characterized 
by a the smearing strength e, smearing level Ns 
and loop size Nl [7]. 

We also performed a calculation using Coulomb 
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gauge smeared operators. We define smeared link 
variables f7/(a;) for s > 1 by the average 

^^W = 77TIP ^ U.{x+pj + qk), (1) 

^ ^ 0<p,q<s 

where j and k are the two space directions or- 
thogonal to i. Here the smeared U^{x) is not 
projected to SU(3), as is done in the case of the 
gauge invariant operators. The product of s se- 
quential U^{y) 

V:{x) = U!{x)...U![x + {8-1)1] (2) 

is then used to construct s x s loops. Therefore 
each operator is uniquely specified by the para- 
meter s. 

At larger (3 we find the gauge invariant oper- 
ators couple more efficiently to glueball ground 
states than do the Coulomb gauge operators. 
Even at smaller /3, however, the required gauge 
fixing makes Coulomb gauge operators computa- 
tionally more expensive. 

Two vacuum subtracted glueball operators Os 
and Os', with the same quantum numbers J-^*^ 
and different smearing parameters r and r', gen- 
erate a propagator Crr'{t) that at large time sep- 
arations approaches the asymptotic form 

ZrZ*,{exp{-mt) + exp[-m{T - t)]}. (3) 

Here T is the lattice temporal period, m is the 
mass of the lightest state, and Z^ is the projec- 
tion {p,\Or\J^'^)- For Coulomb gauge smearing, 
r is the size parameter s, while for gauge invari- 
ant smearing r represents the triple {e, Ns, Nl). 
The Coulomb gauge data was analyzed by fitting 
the propagators obtained from different opera- 
tors to Eq. (3) all at once over a range of sev- 
eral time slices. All the fits were performed using 



Table 1 

Summary of our calculations and results. 



/3 


Lattice 


Count 


Skip 


amo++ 


am2++ 


5.70 


16^x24 


6,050 


400+ (1) 


0.964(42) 


— 


8,094 


10+40 (8) 


0.983(40) 


— 


5.83 


20^x30 


4,002 


400+ (1) 


0.858(43) 


— 


5.93 


24^x36 


4,004 


400+ (1) 


0.811(33) 


1.144(107) 


16^x24 


30,640 


5+20 (16) 


0.786(12) 


1.266(36) 


12^x24 


48,278 


5+20 (16) 


0.771(10) 


1.202(26) 


6.17 


322 X 30x40 


2,005 


400+ (1) 


0.489(31) 


0.816(119) 


24^x36 


31,150 


5+20 (16) 


0.582(10) 


0.828(26) 


6.40 


322 X 30x40 


2,002 


400+ (1) 


0.415(43) 


0.504(61) 


25,440 


5+20 (16) 


0.433(11) 


0.636(23) 



the full correlation matrix for statistical errors. 
For the gauge invariant data we fitted only diag- 
onal propagartors Crr{t) to Eq. (3). The most 
probable combined mass was then found using 
the statistical correlation matrix among the fit- 
ted masses for different r. The required correla- 
tion matrix, and the statistical errors in all fits, 
were found by the bootstrap method. 

Table 1 summarizes our calculations and re- 
sults. For each (3 the first row refers to Coulomb 
gauge operators and the second, if present, to 
gauge invariant operators. The smearing strength 
parameter e was set to 0.25 at /3 = 5.7, and 1.0 
in all other cases. The table also lists the num- 
ber of measurements and the number of 5{7(3) 
Monte Carlo updates between each measurement 
in the format heat bath sweeps + microcanon- 
ical sweeps [8]. Successive measurements in all 
cases appeared to be nearly independent statisti- 
cally. In some cases we found a slight increase in 
error bars if successive propagators were binned 
together to produce a ensemble with fewer mem- 
bers and smaller correlations. In parentheses we 
indicate the bin sizes, generally taken much larger 
than needed, used in our final evaluation of error 
bars. The two different smearing methods gener- 
ally produced statistically consistent mass results. 
We found, however, some discrepancies at the two 
largest /3, presumably due to the small Coulomb 
gauge ensembles yielding unrealiable estimates of 
the size of their errors. Our masses from the 
largest ensembles are all consistent with recent 
calculation by other groups at similar values of (3 



[9,10]. Figs. 1 and 2 show typical effective mass 
plots. A plateau is quite clear in the 0"*""*" data, 
and present, but somewhat more ambiguous, for 

the 2++. 




Figure 1. Effective mass plot for the 0"*""*" glueball 
at /3 = 6.4 for the operator with smearing para- 
meters (e, Ns, Nl) = (1.0, 8, 9). We also show the 
fitted mass and the fitting range. 

The distance of our masses from the infinite 
volume limit can be estimated using Luscher's 
formula [11] 



, . , .r^ ea;p(-V3z/2) 
myz) = to(oo){1 — g 1, 



(4) 



where z = mL and L is the lattice period. From 
our data at /3 = 5.93, we found ^0++ = 682 ± 434 



and 52++ = (1.72±0.85) x 10^, consistent with the 
estimate of Ref. [12]. The gauge invariant data at 
the three largest /3, which we used for the contin- 
uum limit extrapolation, all have Zq++ > 12 and 
Z2++ > 19, which imply masses within 0.5 % of 
their infinite volume limits. This error is signifi- 
cantly smaller than the statistical uncertainty in 
each mass. 
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Figure 2. Same as Fig. 1, but for the 2+"'" glueball. 



To set the lattice scale a we used A^— — from 

MS 

Ref. [1]. Both TOpO and m^a were found to obey 
asymptotic scaling within statistical errors above 
(3 of 5.7. The continuum limit of (^hX-L-a) / {m^a) 



i(o) 



gives Airi- = 243.7 ± 6.8 MeV, consistent with 

° MS 

independent estimates [13,14]. Since the gauge 
part of the QCD action has 0{a^) corrections, 

we extrapolated mass ratios linearly in I ah!^-L- J , 



MS) 

and obtained the continuum limits TOo++/Ai-^ = 

7.14 ± 0.21 and m2++/A^ = 9.68 ± 0.45. Ex- 
trapolations are shown in Fig. 3. Combining the 
errors in ratios and in A^-^, we obtained TOn++ = 

MS^ " 

1740 ± 71 MeV and m2++ = 2359 ± 128 MeV. The 
mass of the observed /o(1710) is close to our esti- 
mate, making it our preferred glueball candidate. 
The 2+"'" is close to at least two well established 
resonances, /2(2300) and /2(2340), which can not 
be distinguished within our statistical uncertain- 
ties. 

We would like to thank Frank Butler for writing 
some of the analysis software which we used, and 
Mike Cassera, Molly Elliott, Dave George, Chi 
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Figure 3. The continuum limit extrapolation for 
both the 0++ and 2++ glueball. 
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